Control of mixed-state quantum systems by a train of short pulses 
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A density matrix approach is developped for the control of a mixed-state quantum system using 
a time-dependent external field such as a train of pulses. This leads to the definition of a target 
density matrix constructed in a reduced Hilbert space as a specific combination of the eigenvectors 
of a given observable through weighting factors related with the initial statistics of the system. A 
train of pulses is considered as a possible strategy to reach this target. An illustration is given by 
considering the laser control of molecular alignment / orientation in thermal equilibrium. 



PACS numbers: 33.80.-b, 32.80.Lg, 42.50.Hz 
I. INTRODUCTION 

The control of molecular dynamics by laser pulses 
has been a long-standing goal in the field of molecular 
physics, both from experimental and theoretical points of 
view (for a recent review, see and references therein) . 
Molecular alignment or orientation induced by an in- 
tense laser field is one of the most challenging processes, 
with applications extending from chemical reactivity to 
nanoscale design 0, 0, 0, EL IS El in this context, 
many efforts, based on physical intuition or on optimal 
control strategies, have been made to construct control 
schemes, but with a limited success. For this reason, 
several works 0, H H E3 have pointed out the connec- 
tion between these results and tools of control theory. In 
particular, they have established kinematical bounds and 
determined the controllability of both pure and mixed- 
state molecular systems. More precisely, since the time- 
evolution operator is a unitary operator, only target 
states (corresponding to a pure or a mixed-state) which 
are unitarily equivalent to the initial state are kinemati- 
cally attainable. This constraint is crucial since it deter- 
mines the upper and lower bounds of the time-evolution 
of the expectation value of an arbitrary observable. How- 
ever, it is important to remark that kinematical attain- 
ability is not synonymous with dynamical realizability, 
i.e. the possibility of finding a path of unitary operators, 
which satisfy the time-dependent Schrodinger equation 
and bring the initial state to the target state within a 
class of physically realizable interactions. For control- 
linear processes (i.e. for Hamiltonians depending in a 
linear way on the control functions) of non-dissipative 
finite-level quantum systems, the dynamical realizability 
of the process depends on the dimension of the Lie alge- 
bra generated by the unperturbed Hamiltonian and the 
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coupling operator 0, • 

We have recently proposed a laser control strategy 
jllL IT^ | , for a molecular system initially in a pure quan- 
tum state, aiming at the maximalization (or minimal- 
ization) of the expectation value of a given observable, 
which does not commute with the field-free Hamiltonian. 
This control scheme consists in applying sudden pulses 
each time a given quantity (such as the expectation value 
of the observable under consideration) reaches its maxi- 
mum. We notice that there are several other methods for 
controlling molecular dynamics using, for instance, adi- 
abatic passage techniques Q, ^3 , factorizations of 
unitary operators or optimal control schemes 0, 0] • 
The main purpose of the present paper is to extend our 
control scheme to mixed-state quantum systems. Our 
concern is twofold : the identification of a well defined 
target density operator and the determination of a con- 
trol scheme for reaching this target. Both steps are very 
general in their principle and their usefulness, such that 
the overall procedure may be extended to different con- 
trol objectives. Alignment / orientation dynamics of a 
diatomic molecule are taken as examples, emp hasi zing 
the effects of non-zero temperature Il9l l2pt |21( . 
The drastic decrease of alignment / orientation with tem- 
perature basically results from an initial superposition of 
rotational states with m ^ (to being the projection of 
the total molecular angular momentum on the field po- 
larization axis), which tends to misalign / misorient the 
molecule. 

The identification of a target state as presented in this 
paper, together with the control strategy to reach it, is 
actually a very powerful tool. The analysis of the target 
state by itself may be also very instructive. An inter- 
esting counterintuitive example is the alignment / orien- 
tation of a diatomic molecule using differently polarized 
lasers. The simple examination of the two target states 
attainable by linear (no modification of to) or elliptical 
(inducing modifications of to) polarized lasers shows that 
they lead to alignment / orientation properties that are 
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very close to each other, rendering the use of elliptical po- 
larizability unnecessary for diatomic systems. As for the 
control strategy, for non-dissipative systems (considered 
in this paper) the perturbations, described by unitary 
operators applied at given time intervals, yield spectacu- 
lar results : the target state may be quite accurately and 
robustly reached within only a few perturbations. 

The paper is organized as follows : we outline the prin- 
ciples of the control strategy in Sec. |n] Section IIIII is 
devoted to the control of the alignment / orientation pro- 
cesses, taken on a parallel footing. Concluding remarks 
and prospective views are presented in Sec. IIVI Some 
details of the discussion are reported in Appendixes lAl 
and El 



II. PRINCIPLES OF THE CONTROL SCHEME 

This section outlines the principles of our control 
scheme in terms of a target state and the strategy to 
reach it. This has roughly an overall similitude with the 
scheme followed for a pure quantum state HJ[ll, but 
the two approaches have to be contrasted at the level 
of the target state. In the pure-state case, it is a spe- 
cific linear combination of rotational states, constituting 
a wavepacket evolving through the Schrodinger equation. 
In the mixed-state case, the target is a density opera- 
tor whose evolution is determined by the Von Neumann 
equation. The specific properties of the target density 
matrix cannot be inferred from a Boltzmann superposi- 
tion of wavefunctions, separately time-propagated, i.e. it 
cannot be obtained from the analysis of Refs [ill fl2| . 

More specifically, we consider a mixed-state quantum 
system whose elements, the density operators p(t), act on 
the Hilbert space TL. A density operator can be written 
as : 



(i) 



where the w k 's are the eigenvalues of p which fulfill the 
following conditions : < w k < 1 and "^ k w k = 1. The 
|^fc)'s are the normalized eigenstates of p. Using atomic 
units (H = 1), pit) evolves according to the Von Neumann 
equation : 



-p(t) = i[p(t),H(t)} 



(2) 



where H(t) is the Hamiltonian of the system. The con- 
trol being exerted through the application of a time- 
dependent external field, H(t) is taken of the form 



H{t) = H + v(t)Hi 



(3) 



where Hq and Hj are, respectively, the field-free Hamil- 
tonian and the coupling operator. v(t) is a real control 
function that can be turned on or off at will. For a given 
self-adjoint operator 0, which does not commute with 



Hq, and an initial state po, the goal of the control is to 
maximize or minimize the expectation value 



(0(t)) = Tr[Op(t)} 



(4) 



More precisely, the goal is to find a path of unitary oper- 
ators U(t), which satisfy the time-dependent Schrodinger 
equation [with an Hamiltonian of the form and bring 
po to the target state pf, for which (O) is a maximum 
or a minimum. 

The set up of the control strategy consists of four steps. 

The first step is the reduction of the original phys- 
ical Hilbert space TC (possibly infinite) to a finite N- 
dimensional subspace TL^ N \ From the physical point of 
view, this reduction can be qualitatively justified by the 
fact that moderate perturbations, i.e. a finite number of 
applied pulses with moderate amplitude, can only trans- 
fer finite amounts of energy to the system, which stays 
thus essentially confined in a finite-dimensional subspace 
[Tl| . This reduction drastically simplifies the study of 
both kinematical constraints and dynamical realizability 
as we will show below. In the subspace THS N ' , we consider 
the reduced operator 0^ N \ which is defined by 



(N) = p{N) 0p (N) 



(5) 



pW being the projector on TL^ N \ To simplify the nota- 
tion, we omit the superscript (N) for the operators Hq 
and Hj in the rest of the paper. Finally, we also as- 
sume that the field-free evolution is periodic with period 
Tq, which means that all the fundamental frequencies of 
Ho in are commensurate and integer multiples of a 

fundamental frequency lu = 2tt/Tq. 

We next define the initial state. We assume that the 
system is initially in thermal equilibrium with a bath at 
temperature T. The initial density operator po is then 
the canonical density operator, written here in the basis 
of the eigenvectors |n) of Hq : 



1 



Po 



■ exp 



Hp 



1 N 

-Y 

Z ^ 



exp 



\n)(n\ , (6) 



where Z = J2 n =i ex P["fcJr"] * s ^ e P ar thion function, fcs 
the Boltzmann constant, and Hq = J2n=i E n \n)(n\. The 
difference of treatment between pure and mixed-state 
quantum systems at this first step is that, while the ini- 
tial state of a pure quantum system is usually a single 
eigenstate of H 0, that of a mixed-state system is a 
density operator involving a thermal average of all eigen- 
states of Hq- 

The second step consists in determining a target state 
Pf which both maximizes (or minimizes) the expecta- 
tion value and is unitarily equivalent to the initial 
state. The general characterization of this target state is 
adapted from a recent work of Girardeau and co-workers 
0, |j • It can be shown that the density operator which 
maximizes (OW) commutes with £>W [jj. These two 
operators can therefore be simultaneously diagonalized. 
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Moreover, due to the constraint of unitary evolution, po 
and pf have the same eigenvalues lo^ with the same mul- 
tiplicity. The expectation value (O^) = Tr[£)Wp F ] 
can thus be written as follows : 

N 

Tr[OWp F ]=J2x k u>a k , (7) 
fe=i 

where the x^'s are the eigenvalues of and the ui ak 's 
the optimal permutation of the eigenvalues lo^s of po, 
which is defined as the permutation which associates the 
highest eigenvalue of pp with the highest one of 0^ N \ 
the second highest with the second one, and so forth. We 
denote pp^ the resulting density operator which maxi- 
mizes and which is precisely our target in the re- 
duced space H.( N \ This density operator is not unique 
if the spectrum of 0^ N ' is degenerate. Here again, we 
emphasize that the target to be reached in our control 
scheme of mixed-state quantum systems is a specifically 
build density operator, as opposite to the case of a pure- 
state quantum system where the target corresponds to 
the ei genv ector of O'"' associated with the highest eigen- 
value |12| . 

The third step deals with the dynamical realizability 
of the target state. In other words, an initial density op- 
erator and a target being given, the question is to know 
if there exists a path of unitary operators that bring the 
system from one to the other. A sufficient condition to 
ensure the dynamical attainability of the target is the 
complete controllability of the system. We recall that a 
control-linear system is completely controllable if the dy- 
namical Lie algebra C generated by the skew-hermitian 
operators Wq and iH\ is u(N). However, it is noted that 
this condition is a strong requirement which is not neces- 
sary. A less strong condition is the notion of connectiv- 
ity where a non-zero amplitude of the evolution operator 
between the initial and the target states only is required 
|22|. Here, in order to clarify the discussion in the gen- 
eral context we assume the complete controllability of the 
process. 

The last step of the control strategy consists in ex- 
plicitely finding a control scheme which allows us to at- 
tain this target. We devise two control strategies which 
are in th e sp irit of our analysis for pure-state quantum 
systems Il2l |23| , while they differ in the definition 
of the targets. More precisely, they consist respectively 
in applying pulses each time {0^ N \t)) (strategy 51) 
or (p [ F N) \p(t)) = Tr[^%(t)]/Tr[(pW_) 2 ] (strategy 52) 
reaches a global maximum within a period To of the field- 
free dynamics. The unitary operator U(t) satisfying the 
Schrodinger equation can then be written as a product 
of terms of the form : 

exp[-iH At]U A , (8) 

the first factor exp[— iHoAt] corresponding to the field- 
free motion and the second one to the application of the 
pulse. At is here the time delay between pulses and A a 



real control parameter. Moreover, we assume that U A 

commutes with £)W or with ppf (depending on the 
strategy 51 or 52) so that its application does not al- 
ter (OW(t)) or (pp\p(t)). More precisely, denoting the 

operator B as O' 9 ' or p^, we have 

Tv[ P (t)B] = Tr[U A p(t)UfB] , (9) 

which proves the claim. As opposed to the value of the 
function Tr[p(t)B], its derivative can be altered by the 
application of the pulse. Indeed, the slope given, before 
the excitation, by 

j t Tr[p(t)B]\ t0 . = Tv[p(t )[H ,B}} = , (10) 
changes into 

j t Tr[p(t)B}\ t0+0 = Tt[p{t Q )Ul l [H Q ,B]U A ] , (11) 

when a pulse is applied at time to- 

We next notice that the two functions (©'"'(t)) and 

(p^p\p{t)) are periodic, continuous and differentiable 
(under field- free evolution) , leading thus to an increasing 
(possibly constant), bounded and therefore convergent 
sequence of maxima. The limit of this sequence is a fixed 
point, i.e. a density operator pj such that 

Tr[p/S] = Tr[U A p f UfB] , (12) 

is a global maximum within a period To, for any value 
of A. In other words, for the mixed state pf this means 
that one cannot find a value of A such that the dynamics 
passes through (under free evolution) a maximum strictly 
larger than the one just before the pulse. We denote by 
T the set of the fixed points. Due to the difficulty to 
determine T exactly, we consider a larger set 5, defined 
as the set of density operators p which fulfill the following 
requirements : 

Tr[p[H ,B}] = , (13) 

and 

Tr[pU^[H o ,B}U A }=0 , (14) 

for all values of A. If the slope of Tr[p(to)B] at time to 
undergoes a change from zero before the pulse : 

^Tr[p(i)S]| to _ = Tr[p(i )[#o,£]] = , (15) 

to a finite value : 

jTT[p(t)B]\ to+0 = Tr[p(to)Uf[H ,B}U A } ± , (16) 

then, Tr[p(t)0] being periodic and continuous, it will 
reach a maximum that is strictly larger within a pe- 
riod Tq- In this way, we can therefore conclude that 
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p{to) £ S and that S contains T. Straigthforward calcu- 
lations (see appendix lAl for details) show that a density 
operator p S S if p commutes with B. We can then ask 
the conditions on Hq, Us and Q( N > such that these den- 
sity operators be the only elements of S and the unique 
fixed points of the strategy. Appendix provides a suffi- 
cient condition on these operators. However, this condi- 
tion does not completely justify the control strategy and 
several open questions are in order : what is the domain 
of attraction of each fixed point, how large is the domain 
of the target state.... 

As a conclusion, we claim that the formulation of this 
control scheme for mixed-state quantum systems is gen- 
eral and powerful in the sense that it emphasizes the 
mathematical definition of a target density operator that 
maximizes the expectation value of the observable O and 
a systematic strategy to reach it through a set of unitary 
perturbations applied at determined time intervals. In 
particular, the construction of the target constitutes the 
real originality of the previous scheme in comparison with 
the one for pure-state quantum systems. The application 
of this control scheme to molecular alignment / orienta- 
tion dynamics will be the subject of the next section, 
where the definition of the target will be highlighted and 
the efficiency of the preceding strategies will be shown 
numerically. 

III. CONTROL OF THE MOLECULAR 
ALIGNMENT / ORIENTATION DYNAMICS 

This section is devoted to the application of the strat- 
egy for controlling the alignment / orientation dynamics 
of polar diatomic molecules. 

A. Description of the model 

We consider a molecule described in a rigid-rotor ap- 
proximation interacting with a linearly polarized electro- 
magnetic pulse. The following molecular Hamiltonian 
H = eJ 2 — V a - (s), already investigated in a previous 
paper [T^ dealing with the control of pure-state quan- 
tum systems, will serve as an illustration. The time- 
dependent Schrodinger equation which governs the dy- 
namics can be written as : 

*§-/(0, cf>; s) = [eJ 2 - V a , (sM(e, cf>; s) , (17) 

where J 2 is the angular momentum operator, 9 and (ft 
denoting respectively the polar angle and the azimuthal 
angle. The operators V a;o , corresponding to the radia- 
tive interaction terms for the alignment / orientation pro- 
cesses, are defined respectively as follows : 

V a (s) = E 2 (s) cos 2 9 + F 2 (s) , (18) 

and 

V (s) =E (s) cos 9 , (19) 



E a - a and F a being the dimensionless interaction strength. 
These time-dependent functions are equal to : 

El(s) = Aarf 2 (rs)/2 (20) 
F!(s) = a ± rf 2 (rs)/2 (21) 
E (s) = Mor/(r S ) , (22) 

where fj,o is the permanent dipole moment, Aa = a\\ —u±_ 
is the difference between the parallel ot\\ and perpendicu- 
lar aj_ components of the polarizability tensor and f(t) is 
the enveloppe of the laser pulse. Note that Eq. I|17(l is ex- 
pressed in terms of the dimensionless parameter e = tB 
and s — t/r, where B and r are, respectively, the ro- 
tational constant and the pulse duration. The reader is 
referred to Ref. [12j for the explicit derivation of Eq. 
(JT7J. We also recall that, due to the cylindrical sym- 
metry, the projection m of the total angular momentum 
j on the field polarization axis is a conserved quantum 
number. 

From the practical point of view, we are interested in 
short duration pulses, the duration of the pulse being 
compared with the molecular rotational period T rot = 
ir/B. Such kicks allow us to consider a sudden approxi- 
mation and to derive a simple effective evolution opera- 
tor describing the interaction with the pulse pil |25| . We 
have 

U a = exp[iA a cos 2 9} , (23) 

where A a = E 2 {s)ds (the contribution of F 2 (s) has 
not been considered as it corresponds to a pure phase 
factor) for the alignment process and 

U = exp[iA cos6>] , (24) 

with A a = J Q E a (s)ds for the orientation. Using these 

latter propagators and the fact that [Oi%' , p^D\ — 0, one 
can easily check the necessary commutation relations 

[U a;o ,oW} = lU a . iO ,pi$}=0, (25) 

O a -o being defined as O a — cos 2 9 and = cos 9. 

The efficiency of the alignment / orientation is charac- 
terized by a thermal average of O a;o over the rotational 
levels : 

(O a;o )(s)=Tr[p(s)O a;o ], (26) 
with as an initial condition 

P° = \ E E \j,m)e- BjU+1)/kBT {j,m\ , (27) 

m6Z j>|m| 

where Z = J2 m ezY,j> m e ~ Bj( ' j+1)/kBT is the partition 
function and the \j, m)'s the eigenvectors of J 2 . It is to be 
noted that the goal of the control of molecular alignment 
/ orientation is not only to obtain large absolute values 
of (O a ;o) but also a persistance of this effect. We will 
see that one can increase the duration of the alignment / 
orientation by choosing adequately the dimension of the 
molecular rotational space while keeping a high efficiency 
of the process. Details on this point will be given in Sec. 

mm 
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B. Kinematical bounds and dynamical realizability 

We first analyze the kinematical bounds of the pro- 
cess. For doing so, we have to reduce the physical Hilbcrt 
space Ti to a finite subspace H^ maa \ jmax being the high- 
est j's for which the corresponding rotational levels are 
significantly populated. Simple algebra shows that the 
dimension N of ft tf j s 



N = {jmax + l) 2 



(28) 



We recall that the density operators p^"^ associated 
with the kinematical bounds are given by [Eq. Q] : 



Jim. 
Popt 



N 

E 



^ k \x {k) )(x (k) \ 



(29) 



the Wk's being the eigenvalues of the statistical density 
operator [Eq. 1(771) ]. a a permutation of these eigenvalues 
depending on the set of density operators which are dy- 
namically attainable and \x^) the eigenvectors of oif?. 
If we assume the complete controllability of the system, 
is a maximum for the permutation which asso- 
ciates the highest eigenvalue of p^' with the one 

(j ) 

and so forth. We call this maximum, Pop7 the optimal 

(j ) 

density operator. It is noted that Pop7 is not unique 

as the eigenvalues of O^J are degenerate. 
However, a diatomic molecule driven by a linear polarized 
laser pulse is a system which is not completely control- 
lable, to being a good quantum number, the system is 
decoupled and the Hilbert space ft^™^) can be written 
as the direct sum of the dynamically invariant subspaces 
ft 



Un 



^ u m——j n 



(30) 



Moreover, as the subspaces of Tim ma ^ of a given parity 
of j are not coupled by the operator O a N ^ , another sub- 
division has to be considered for the alignment process 
for which we can write 



ftCimaaO 



m.ev. 



ft 



(jm« ) ^ 

od. I 



(31) 



The subscripts ev. and od. correspond respectively to the 
even and odd values of j. To simplify the notation, wc 
omit these subscripts below when confusion is unlikely. 

Let P m be the projector onto the subspace Tim ma ^ 
and the observables pffl and 0^ be the restrictions 



c \ We denote by wjj, 7 ™' and 



of P W and OT } on ft}„ - > — ■ 

(m) (j ) 

\ k the eigenvalues of these two operators in ft™ a . 
As all the subspaces are initially populated, the optimal 
density operator for a linear polarization p\{™ a ^ is the 
direct sum of the optimal density operators [with this 
time the general definition given by Eq. Q] p^ in each 



subspace ft 
follows 



^[O^pt 



Tv[0 ( 



W Jim. 

Plin 



] can be written as 



E E 



(m) (m) 
J k Xfc 



(32) 



fe=l 



where the 's and the xi™"* ' s are ordered for each value 
of m. Note that a similar result taking into account the 
parity of j can be established for the alignment. 

We now examine more closely the different aspects 
of the controllability of these processes. We begin by 
analysing the individual controllability of each subsys- 
tem. Within each decoupled invariant subspaces corre- 
sponding to different values of m, the controllability of 
the subsystem is related to the dimension of the dynami- 
cal Lie algebra generated by the free Hamiltonian and the 
coupling operator. We consider a type of Hamiltonian 
which has already been investigated by S. G. Schirmer 
and co-workers [jjj, |26( and it was shown that each sub- 
system is completely controllable together with the re- 
striction on the parity of j for the alignment. 

The density operator p\ ™ a * is dynamically realizable 

if all the subsystems acting on Jim are comple tely 
controllable, but this condition is not sufficient |27l |28| 
to ensure the controllability of the whole system. In- 
deed, it has not been proved that all the subsystems can 
be simultaneously controlled by the same control pulse. 
Appendix [5] provides numerical and theoretical results 
that indicate that the whole system is not simultaneously 
controllable. However, a less strong notion of simultane- 
ous controllability can be derived by taking into account 
the particular symmetry of the system, i.e. the fact that 



its dynamics is the same in the subspaces Ti. 



Umax) 



and 



n 



Finally, with this latter condition on the dy- 



namical Lie algebra (see Appendix IB1 for details), it can 
be shown that the target states are dynamically attain- 
able for the alignment / orientation processes. 

Returning to the description of the kinematical 
bounds, Figs. ^ an( i [21 display the maximum efficiency 
for the alignment / orientation processes of the molecule 
LiCl as a function of j m ax- A very striking observa- 
tion is the amount of alignment / orientation that can 
be achieved by the strategy using solely linear polarized 
pulses, namely within 10% of the optimal kinematical 
limit. Such an agreement is expected in the high rota- 
tional excitation limit where sin# = m/j tends to zero, 
producing alignment even though the initial set of quan- 
tum numbers m has not been modified by the radiative 
field. Which is more unexpected is that such a circum- 
stance occurs for low excitation (i.e. jmax — 8). This 
leads to the important conclusion that alignment / ori- 
entation control could be satisfactorily achieved referring 
to linearly polarized pulses. This, as opposite to more 
sophisticated schemes working, for instance, with ellipti- 
cally polarized pulses, may present an interest from an 
experimental point of view. We also notice that the op- 
timal alignment / orientation increases as jmax increases 



max 



max 



FIG. 1: Maximum alignment efficiency as a function of jmax 
(see text) for the molecule LiCl in the case T — 10 K (crosses) 
and T = 5K (open circles). Solid and dashed lines, which are 
just to guide the eye, correspond respectively to the optimal 
maximum and the one for a linear polarization. 




J max 

FIG. 2: Same as FigQ but for orientation. 




max 



FIG. 3: Maximal alignment efficiencies that can be obtained 
with a linear polarization (crosses) and associated duration 
(open circles) as a function of jmax (see text) for the molecule 
LiCl in the case T = 5K (dashed lines) and T = 10K (solid 
lines). 



FIG. 4: Same as Fig [2] but for orientation. 



or the temperature T decreases. However, it is impor- 
tant to realize that a larger N corresponds to a better 
efficiency, but leads to a shorter duration. This point is 
shown in Figs. |3and0]where At/T rot is the relative dura- 
tion of the field-free alignment / orientation over which 

^[Plin 'olio] remains larger than 0.5. We then see 
that a compromise has to be made between maximum 
efficiency and duration [Tl| . For the alignment / orien- 
tation of the molecule LiCl, we choose j max — 8 which 
allows us to obtain a duration of the order of 1/20 of 
the rotational period and an optimal efficiency of 0.8 at 
least for T = 5 K. It is finally very important to note 
that only poor alignment / orientation efficiency and du- 
ration can be achieved for high temperatures whatever 
the sophistication of the excitation scheme be (different 
polarization directions, for instance). About 25 to 30% 
of the efficiency is lost when heating from 5 K to 10 K. 



C. Control of the dynamics 

The results are presented for the molecule LiCl with 
the rotational temperature T = 5 K. As these results 
only depend on two dimensionless parameters t/T rot 
and B/(ksT), they are transferable to any particular 
molecule. Numerical parameters are taken to be : A a = 
2, A a — 2 and e = 0.01. This amounts to a pulse duration 
of about 0.2 ps and a field amplitude of 1.5 x 10 5 V.cm -1 
|24| for the case of orientation. The rotational constant 
B = 0.70652 cm -1 (value at the Li-Cl equilibrium dis- 
tance) is chosen so as to reproduce the principal features 
of the diatomic molecule LiCl. In the different figures, 
time is indicated in fractions of the LiCl rotational pe- 
riod. 

In this section, we are interested in the dynam- 
ical behavior of the two quantities Tr[(9 a;Q p(t)] and 

Tr[p[ 4 8 ? |p(t)]/Tr[(p[ 4 8 ? |) 2 ], which correspond respectively to 
the expectation value of the observable O a -o and to the 
projection on the target state p,?|- Except for Figs. |H1 it 
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FIG. 5: Alignment dynamics during and after the train of 
pulses determined by strategy SI, with pulses acting at the 
global maxima of Tr[p(t) cos 2 8] : panel (a) for Tr[p(i) cos 2 8] 
and panel (b) for Tr[pj2p(t)]. The solid line corresponds to 
the exactly propagated density operator p(t) and the dashed 
line to the propagation of pit) in the subspace . The train 
of pulses is displayed on panel (b), the optimal alignment 
and the one for a linear polarization in Tt^ are, respectively, 
indicated by the horizontal solid and dashed lines on panel 
(a). 



is noted that all results are displayed for times starting 
at the first pulse (taken thus at the origin) and extending 
up to one rotational period T rot after the last pulse. This 
helps in showing the complete field-free behavior of the 
dynamics and the well-known revival structures. 

The results of the strategy 51 for the alignment and 
orientation processes are illustrated in Figs. [5] and [S] 
As could be expected, the alignment dynamics presents 
a more strongly oscillating structure as compared with 
the orientation one. This difference explains the differ- 
ent choice of parameters A, N and the number of kicks in 
these two processes. In all cases, panels (b) show that the 
density operator p(t) gets close to the target state p®, 
which also shows that the rotational dynamics roughly 
resides within TiS 9 ^ . It can furthermore be checked that 
the postpulse dynamics lead to very good results in terms 
of efficiency and duration. Indeed, an efficiency of about 
0.75 and a duration of about 1/20 of the rotational pe- 




FIG. 6: Same as Fig. but for the orientation dynamics. 

riod (that is about 1 ps for a light molecule like LiCl) is 
obtained for the orientation. To our knowledge, this is 
the largest duration and efficiency achieved up to date 
for such a thermal ensemble. 

In Figs. [7| we display the mean values Tr[0 a;o p(i)] 

and Tr[/9^/3(t)]/Tr[(p^) 2 ] computed with the strategy 
52 for the case of orientation. It is seen that very close 
results are obtained with respect to Figs. [fJland the strat- 
egy 51, which shows the similarity between the two con- 
trol schemes. Note that a smaller number of pulses, i.e. 
9, is used for 52. 

Finally, the evolution after the last kick for the align- 
ment / orientation processes are plotted in Figs. |H1 
which, in addition, also display the dynamical behavior 

of the two optimals, p u ^ and p K Q J t . We notice that the 
control strategy allows us to obtain a density operator 
with free dynamics very close to the optimal one, both in 
efficiency and duration. Once again, these results stress 
the importance of the target states since they give a clear 
insight into the free dynamics of the molecule. 



IV. CONCLUSION 

A precise mathematical definition of a target provides 
a completely new insight to the control issue of a system 
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FIG. 7: Same as Fig. El but referring to the strategy S2, 
with pulses acting at the maxima of Tr[p^/j(t)]. 



using a time dependent external field. Complementary to 
what is done for pure-state quantum systems |12| where 
the target is an eigenvector of a given observable pro- 
jected onto a finite dimensional Hilbert subspace, the tar- 
get in the mixed-state case is a non-trivial density matrix 
constructed as a specific combination of the eigenvectors 
of a given observable through weighting factors related 
with the initial statistics of the mixed system. Once the 
target is identified and its attainability established, two 
possible strategies to reach it are thoroughly analyzed. 
They amount to apply a series of unitary perturbations 
each time Tr[p(t)0^ N '] or Tr[p(t)p opt ] reaches a global 
maximum (or minimum) under free evolution. 

An illustration is given by considering the laser control 
of molecular alignment / orientation in thermal equilib- 
rium. Although some previous works have already used 
the strategy of applying trains of short laser pulses, the 
originality resides here in the practical construction of the 
target density matrix and in comparing the optimal align- 
ment / orientation actually obtained to the ones which 
are the best theoretically possible within the reduced sub- 
space supporting the dynamics. 

As a prospect, an important open question in this field 
is the applicability of the present method to more compli- 
cated systems involving, for instance, an interaction with 



FIG. 8: Postpulse alignment [panel (a)] and orientation 
[panel (b)] dynamics of the molecule LiCl after interaction 
with a train of short pulses (the time t = corresponds to 
the first kick) . Solid, dashed and dash-dotted lines correspond 
respectively, to the averages calculated with the exact density 
operator, the optimal one and the one obtained with a 

f 81 

linear polarization In these two cases, the strategy SI is 
used. 



a physical environment or a dissipative system (2iJ |3(| . 
The answer is not obvious because the strategy would 
involve non-unitary evolution and special attention has 
to be paid to decoherence 0, Moreover, for an 

important class of dissipative systems such as laser cool- 
ing, the dissipation can increase the purity of the state 
pM IscL Isflj. In this case, th e op timal target state would 
be the one defined in Ref. [Tl| for pure-state quantum 
systems. 
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APPENDIX A: ANALYSIS OF THE SET OF 
FIXED POINTS 

In this appendix, we analyse the set S with the hy- 
pothesis of complete controllability. For that purpose, 
we introduce the vector space V generated by the op- 
erators [H ,B] and U^H^Bp^, where A £ R. One 
can reformulate Eqs. (|T3|) and QTty by stating that p is 
an element of S if p £ V x , where V ± is the orthogonal 
space of V. It is noted that V -1 does not only contain 
density operators. Moreover, it can also be shown that 
p £ S if the density operator p commutes with B. In- 
deed, as [B, ETj] = 0, this remark can be easily checked 
by calculating the following expression : 

Tt[pUf[H ,B}U A ] = Tr[p[UfH U A ,B}} (Al) 
= TT[(Bp-pB)U A 1 H Q U A ](A2) 

At this point, the question that naturally arises is the 
condition on Ho, B and U A such that S only contains 
these latter elements. The condition is the following : 
The density operators which commute with B are the 
only elements of S if and only if the dimension of V is 
N 2 — J2i n h where the n^s (i = 1, • • • , P) are the multi- 
plicities of the eigenvalues of B. 

The proof is straightforward. We note \b i; k)(k — 
1, • • • , m) the eigenvectors of B. A density operator 
such that [p( N \B] — is of the form : 



p 

E 

i=l k.k 



ai;k;k'\bi;k)(h;k' 



(A3) 



where the ai-k-k^s are complex numbers such that the 
eigenvalues cok of p^ fulfill the two conditions : < 
ujk < 1 and Y2k=i UJk = 1- The vec tor space gener- 
ated by these operators is of dimension ^ nf and it 
is a subspace of V x . In addition, [/^ [iJ , 1 be- 
ing a skew-hermitian operator, it can readily be shown 
that C contains V. One finally deduces that the dimen- 
sion of V is at most N 2 (the complete controllability be- 
ing assumed), which completes the proof (we recall that 
dim V + dim V" 1 = N 2 ). 



APPENDIX B: SIMULTANEOUS MIXED-STATE 
CONTROLLABILITY OF THE MOLECULAR 
ALIGNMENT / ORIENTATION DYNAMICS 

In this appendix, we consider the simultaneous control- 
lability of the alignment / orientation processes driven by 
a linear laser field. A sufficient and necessary condition to 
ensure the simultaneous mixed-state controllability of a 
system is given by a condition on the dimension of the dy- 
namical Lie algebra C (generated by the skew-hermitian 

operators i(Ho)( N ^ and iO a N J) which must be equal to 



D = r+(J max + l) 2 -l+2 [{J m ax-m+l) 2 -l\ , (Bl) 



where r is the rank of the matrix 

/ Tr[H ] m =j max Tr[0 O ; O ] m =j mM \ 



T 



Tr[if ] 



Tr[O a;Q ] 



V Tr[ff ] 



Tr[O a;o ] 



(B2) 



/ 



For the case of orientation, straightforward calculations 
lead to the following expression : 



T 



( jmax^jmax H~ l) N 



ESifc(* + i) o 



V jmaxi^jrnax H~ l) J 



(B3) 



which implies that the rank r of this matrix is r = 1. 
In the case of alignment where the trace has to be taken 
separately in the subspace Ti. m ma ^ for the odd and even 
values of j, this rank is equal to r = 2. 

We now analyse the simultaneous controllability of the 
orientation process. We have already mentionned that 
each subsystem for each value of m is completely con- 
trollable. We consider independently the two subsystems 
corresponding to the values m = ±(j max — 1), jmax > 1, 
which can be viewed as two non-interacting two-level sys- 
tems. In this case, we have : 



(0 ) s = da x 9 da x , 



(B4) 



and 



o)s 



diag(4 \E[ >)ediag(£,r J ,£i 



( + ) B-C+h 



(B5) 



where d = + 2j 



and E 



(-) 



E 



(+) 



max ■, 



(-) 



E 



(+) 



0V. 



l)ir, 



{jr. 



l)jr, 



The subscript S in- 



dicates the restriction of the operator to the two subsys- 



tems. AsE[ +) -E ( +) 



{ - ] -E^ 



it can be shown |2f 



that this system and therefore the whole system (when 
jmax > 1) are not simultaneously controllable. Note that 
such a result can also be established for the alignment 
process. 

At this point, it is important to remark that a con- 
straint on the dynamics has not been taken into account 
in the condition of Eq. 1|B1|) . Using the fact that the dy- 
namics in the subspaces 7i m ma ^ and 7i < f™ ax ' 1 is the same, 
a new criterion on the dimension of the dynamical Lie al- 
gebra can be derived for the orientation process, where 



only the subspaces Ti 



with m > are considered : 



D' = r+(J max + l) 2 -l+ V [{Jmax-m+lf . (B6) 
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A similar condition can be established for the alignment 
together with the condition on the parity of j. 

From a numerical point view, we can calculate the di- 
mension of C for low values of j m ax |34| . The results of 
these computations are presented in Tab. HI and ITU 



TABLE I: Dimension of the dynamical Lie algebra £ as a 
function of j max for the orientation dynamics. The third and 
fourth columns indicate respectively the dimensions D and 
D' of £ needed for the general [Eq. 1B1H and restricted [Eq. 
1B6H simultaneous controllabilities of the system. 



TABLE II: Same as Tab. |TJ but for alignment. 



J max 


dim(£) 


D 


D' 


l 


4 


4 


4 


2 


12 


15 


12 


3 


27 


38 


27 



Jmax 


dim(£) 


D 


D' 


l 


2 


5 


2 


2 


5 


16 


5 


3 


11 


39 


11 



The results of Tabs. [I] and [H] suggest the hypothesis 
that the alignment / orientation processes are simultane- 
ously controllable in the restricted sense defined above. 
Moreover, as by construction the target states fulfill the 
constraints of the dynamics, we can conclude that these 
targets are dynamically attainable. 
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